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ABSTRACT 


The specific production of tritium in sampleo of ‘LW and ‘I&I 
embedded 
in a 600-mm-diam sphere of OLiD irradiated by a central source cf 14-MeV 
neutrons has been detmmined by measuring 
the activity of the hydrogen 
evolved from the samples of each isotope at each of five different radii in the 
OLiD assembly. The entire process of decomposing the LiH, transferring 
the 


evolved gas into counters, and determining the decay rate was standardized 
by processing IJH 
samples irradiated 
by thermal neutrons for which the 
%i(n,a) 
cross’ section is well known. 
The specific production of tritium in ‘LiH and TAH (embedded samples) 
and the activation of radiochemical 
detector foils of %c, 
‘Y, 
‘Zr, 
1%, 
1’lIr.~7,‘WIr.8~, lWAU, ‘U, 
and ‘W 
placed at various positions in the ‘LiD 
sphere have been calculated and compared with the experimental data. One- 
and three-dimensional 
Monte Carlo and S. neutron-transport 
calculations 
were 
performed. 
The 
most 
reliable 
(three-dimensional 
Monte 
Carlo) 


calculation is in reasonable agreement with both the tritium-production and 
the radiochemical-activation 
data. 
The 
existing 
discrepancies 
between 
calculation and experiment appear largely 
attributable 
to uncertainties in 
some tritium-production 
and radiochemical-activation 
cross sections. 


.- 
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PART 
I 


MEASUREMENT 
OF TRITIUM 
PRODUCTION 
IN A SPHERE OF 
%iD HUtADIATED 
BY 14-MeV NEUTRONS 


by 


. 


A, Henunendinger, 
C. E. ~an, 
E. R. Shunk, 
A, N. Ellis, and J. M. AUWa 


ABSTRACT 


The opecific production 
of tritium 
in G 606-mmdiam 
sphere of %iD 
ir- 


radiated 
by a central 
source of 14-MeV 
neutrons has been determined 
by 
rnaaouring the activity of the tritium in samples of ‘LiH and ‘LiH contained 
in sealed quartz Gmpulea embedded in the sphere. Results are reported for 
several 
umplea 
of each isotope at each of five different 
radii in the as- 
wmbly. 
The entire process of decomposing 
the LiH sampleu, transferring 
the evolved 
Cas into counters, and determining 
the decay 
rate was stan- 
dardiaad by procemins 
LiH 
samples 
irradiated 
by thermal 
neutrons for 


which the ‘Li(n,a) 
cross section is well known. 


——.————————————————— 


I. INTRODUCTION 


The 
determination of the amount 
of tritium 
produced in a blanket of lithium surrounding a 
neutron source has long been of interest in connec- 
tion with nuclear weapons and controlled ther- 
monuclear reactora. If a sustained thermonuclear 
reaction requirestritium as a fuel, the only plausible 
way to keep such a system operating would be to 
create tritium as a by-product of 14-MeV neutrons 
produced in the reactor. In 1956 Marvin Wyman* 
reported the tritium production in samples of ‘L1 
metal embedded in a 600-mm-diam sphere of LiD 
surrounding a source of 
14-MeV neutrons. His 
results confirmed that tritium is indeed produced in 
substantial quantities by the reaction ‘Li(n,n’a). 
Several measurement 
of tritium production in 
blanketa containing natural lithium are shown to be 
in fair agreement with calculations.s-’ Muir and 
Wyman’ and Qaim, Wtilfle, and St&klin” have dis- 
cussed the need for integral data on tritium produc- 
tion, and they have analyzed earlier experiments. 


We 
report 
here 
an 
experiment, 
similar 
to 
Wyman’s, in which we irradiated an assembly of 
‘LiD 
containing test samples of ‘LiH and 7LiH. 
Teledyne 
Isotopes, 
Inc., 
Westwood, 
NJ 
07675, 
analyzed these samples by evolving hydrogen from 
them and counting the tritium decay by a method 
resembling, but developed independently of, that 
described by Qaim.7 Our experiment, which was 
designed to investigate the transport of neutrons in 
OLiD and tc 
determine the tritium production, 
provides benchmark measurement for comparison 
with calculations. Results for tritium production in 
“ .i in such an assembly have not previously been 
repmted. 


IL 
THE 
‘LiD 
ASSEMBLY 


Nesting hemispherical shells of ‘LiD were con- 
structed to fit around the target assembly of a 
Cockcroft-Walton (C-W) accelerator. The inside and 
outside diameters of the assembly were 44.4 and 600 


. 


. 


. 


2 


mm. Details of dimensions, masses, and isotopic 
composition are given in Table I. Small cavities at 
the shell interfaces held samples of 6LiH and ‘LiH; 
these samples were in sealed quartz ampules (Fig. 1) 
in two sizes: 10-mm o.d. for the two inner-shell inter- 
faces and 18-mm o.d. for the outer shells. The am- 
pules were filled, weighed, and sealed in an at- 
mosphere of helium.* All ampules were taped in 
place in the assembly, as shown in Fig. 2. At least 
two empty ampules (treated in the same manner as 
the filled ampules) were positioned at each radius 
for background measurement. 
The 
neutron-source 
flux 
distribution 
was 
measured over 47rsr using both liquid-scintillator 


*M ‘LID parta and LIH samples were prepared by Union Carbide 
Corporation, Oak Ridge, TN 37830. 


TABLE 
I 


SPECIFICATIONS 
FOR 
‘LID 


HEMISPHERICAL 
SHELLS 


Diameter 
(mm) 
Mass 
%i 


Inside 
Outside 
(kg) 
— 
. 
_ 


44.4 
100.0 
0.1783 


0.151 


102.0 
152.3 
0.449 


0.490 


154.3 
252.0 
2.465 


2.330 


254.0 
400.0 
9.133 


9.430 


402.0 
600.0 
30.000 


29.200 


Total mass = 
83.826 


(at.%) 


95.59 


95.59 


95.59 


95.68 


95.68 


(Wt%) 


94.89 


94.89 


94.89 


95.00 


95.00 


uH 
8 OR16mm 


1-10 OR18mm 4 


Fig. 1. 
Quartz ampule, 
with dimensions 
for two sizes, 
which held the lithium hydride samples for ir- 
radiation. 


Fig, 2. 
Lithium 
hydride 
samples 
taped 
in place 
in 
cavities 
in the 8LiD hemishell. 
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Fig. 3. 
Map of the neutron flux over the sphere around 
the C-W machine 
target. 
The measurements 
made 
wing 
a NE-213 
(Nuclear 
Enterprises, 
Inc.) scintilZator biased at 12 MeVare 
shown as 
circles, and those obtained 
using a cadmium- 
shielded 
‘86(J fission 
detector 
are shown 
as 
triangles. 


and fission detectors. This distribution is displayed 
in Fig, 3 as a function of the azimuthal angle 0 (0° 
being the direction of the deuteron beam) and the 
polar angle O, the angular displacement about the 
axis defined by the deuteron beam (0° being in the 
horizontal plane to the right side looking into the 
beam). The observed flux distribution is not quite 
isotropic because of a superposition of two effects: 
(a) neutron emission from the recoiling compound 
nucleus varies relative to the 90° flux by +5% in the 
forward direction to –57. 
in the backward direction 
and (b) the flux is attenuated by about 10”Ain direc- 
tions in which neutrons traversethe target backing. 


The copper target wafer, which had a surface layer 
of tritium embedded in titanium, was oriented in a 
vertical plane at an angle 45° from the direction of 
the horizontal deuteron beam. The ‘LiH samples in 
the assembly were located in the vertical – 110° 
plane (B-B), the ‘LiH in the vertical –65° plane (A- 
A), as shown in Fig. 4. The location of a sample in 
either plane is specified by its radius r and by a polar 
angle w measured from the parting plane looking 
into the beam. Care was taken in positioning the 
samples to avoid locations at which there might be a 
neutron flux perturbation; thus no two ampules were 
placed along the same radius, near the parting plane 
of the sphere, or in the plane of the C-W target where 
there are flux perturbations as large as +8%. These 
precautions were taken even though variations in 
the flux were probably smoothed out at distances 
greater than one mean free path in the ‘LiD—108 
“mm for 14-MeV neutrons. 
Figure 5 shows the lower half of the assembly in 
place around the target, and Fig. 6 shows the par- 
tially completed ‘LiD assembly. The assembly was 
supported in a 520-mm-diam hole in a 6.4-mm-thick 
circular aluminum plate. The plate was suspended 
from an assembly frame (Fig. 7) constructed to 
mount massive objects around the C-W machine 
target. 


111. STANDARDIZATION 


A. General Procedure 


To 
evolve 
hydrogen 
from 
lithium 
hydride, 
Teledyne crushed each quartz ampule in a stainless 
steel tube that they heated until there was no further 
evolution 
of 
gas. The 
volume 
of 
the 
evolved 
hydrogen was measured, and hydrogen carrier was 
added to the small samples only. The hydrogen gas 
was converted to water, and the water was diluted 
and divided into several parts in measured ratios. To 
assay one of these subdivided samples, they con- 
verted it back to hydrogen and transferred the gas 
into a shielded proportional counter. 
In addition, a portion of the tritiated water from 
each sample was mixed with a liquid scintillator, 
and 
the 
sample 
count 
rate 
was 
determined. 
Teledyne calibrated the counting system routinely 
by counting part of a standard sample, NBS 4925, 


, 


t 


4 


. . 


.. 


Fig. 4. 


Schematic 
showing positions 
of lithium hydride samples in the “LiD assembly. 
cles for some samples 
indicated 
that the stem 
of the ampule 
is toward 
the 
angles above the parting plane are negatiue. 


which they had used for about 
10 years; this 
provided Teledyne’s primary standard. They used 
two additional calibration methods: 
(1) counting a fraction of NBS 4926 standard 
sample, which had been calibrated by the 
National Bureau of Standards on 1July 1976; 
(2) comparing Teledyne’s count on a fraction of 
our sample 327 with NBS’s count on another 
fraction of the same sample. 


The small cir- 
observer. 
The 


A comparison of calibrations 1 and 2 (Table II) 
shows that the Teledyne and NBS determinations of 
tritium activity agree within 1.5Y0,which is within 
Teledyne’s standard deviation. NBS estimates their 
systematic error at 20A;this is the same systematic 
error inherent 
in 
the 
Teledyne 
measurement, 
because both groups used fractions of the same 
tritium standard source. 


5 


\ 


Fig. 5. 
Bottom 
half 
shells 
of 
the 
8LiD 
assembly 
mounted 
around 
the 
target 
of 
the 
C-W 
machine. The tube at the lower left is for target 
insertion and cooling; at upper left is the beam 
entry 
tube, 
and at upper right is the alphu- 
particle 
counter 
tube 
for 
monitoring 
the 
neutron flux. 


Fig. 6. 
Shells 
of 8LiD mounted 
around 
the 
tritium 
target of the C-W machine. 
The bottom half of 
the assembly 
is complete; 
the outer top three 
hemishells 
are missing. 


Because there may be other systematic errors, 
which we believe are associated with gas handling by 
Teledyne (see Sec. I.11.E),tritium production results 
in lithium hydride irradiated by thermal neutrons 
were used to calibrate the whole process of gas 
handling and counting. These lithium hydride sam- 
ples, each containing a gold fluence monitor, were ir- 
radiated in the thermal column of the Los Alamos 
Scientific Laboratory (LASL) Omega West Reactor. 
The ratio of the thermal capture cross section of ‘Li 
to that of I’”Au is well enough known””’ that tritium 
counting by Teledyne of these reactor-irradiated 
samples 
provided 
a calibration 
of 
Teledyne’s 
measurements relative to the absolute counting of 
the gold by LASL. 


. . 


.. 


Fig. 7. 
Completed 
6LiD assembly 
surrounding 
the 
target of the C-W machine. 
The aluminum 
cir- 
cular 
mounting 
plate 
is suspended 
from 
a 
pneumatic 
assembly 
machine 
by three rods. 
Three sample ampules 
can be seen on the out- 
side surface. 


. 


b 
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TABLE 
II 


COMPARISON 
OF TRITIUM 
DETERMINATIONS 
IN WATER 


Sample 


NBS H-3 


LASL-327 


B. 
Backgrounds 


BY 
TELEDYNE 
AND 
NBS 


Teledyne Measurement 
NBS Measurement 
Ref. Date 
[108dis”s-’”(g 
H,O)-’] 
[l& dis”s-’”(g 
H,O)-l] 


8/16/76 
2.242 
2.292 
av=~ 
10/11/76 
1.960 


A surprise came in our investigation of samples of 
LiD, the material that we originally intended to use 
in the ampules, for they had tritium backgrounds 
high compared with the activation expected in the 
C-W irradiations. For completeness, we report these 
background levels in Appendix A. Inasmuch as the 
deuterium abundance in our LiD was typically 99.9 
at.%, it is reasonable to suppose that an increase of 
the deuterium content by a factor of 7000 above the 
natural level will necessarily increase the tritium 
fraction of hydrogen even more; the natural abun- 
dance of tritium augmented by enrichment appears 
sufficient to account for the observed backgrounds. 
In our samples of ‘LiH, the tritium background was 
below Teledyne’s detection limit; the samples of 
7LiH contained appreciably more tritium, but never 
more than a few per cent of the activity induced by 
irradiation. Measured 7LiH backgrounds are listed 
in Table III. 


C. Thermal 
Activations of LIH 


The test samples of lithium hydride irradiated in 
the thermal flux contained a low concentration of CLi 
(about 0.1%) to limit perturbation of the neutron 
flux by absorption in the sample. Even so, the ‘Li 
and ‘H in the samples produced a combined flux 
depression, with approximately equal contributions, 
of about 5Y0, a depression for which we could ade- 
quately correct. 
The gold fluence monitor obviously should have 
been uniformly distributed in the sample, but we 
found no way to do this. Gold powder settled quickly 


2.276 
1.944 


TABLE 
III 


Ratio 
T/NBS 


0.9960 
1.0083 


TRITIUM 
ACIYVITY 
IN UNIRRADIATED 


AMPULES 
OF ‘LiH 


Mass 7LiH 
‘Ihitium Decay Rate 
Sample No. 
(g) 
[dis*s-’”(gLiH)-l] 


125 
127 
128 
179 
176 
182 
189 
190 
316 
322 
323 


0.1467 
0.1589 
0.1694 
1.0125 
1.0069 
0.9501 
1.0028 
1.0027 
0.8982 
1.0683 
1.0871 


2.317 
4.900 
3.250 
2.750 
2.667 
3.260 
3.600 
3.516 
3.583 
8.617 
3.lm 


av = 3.78 + 0.52 


to the bottom of the ampule, and the use of solutions 
promised to add more problems than it solved. In 
each ampule we placed a sample of approximately 
10 mg of 0.05-mm-diam gold wire wound in a helix 
on a 1.6-mm-diam mandrel. The helix was stretched 
axially to span the inside diameter of the ampule 
(about 8 or 16mm for the two sizes of ampules), and 
we attempted ta locate the helix along the axis of the 
fill tube. X-radiographs of the filled ampules showed 
that not all helixes were collinear with the fill tube, 
but we judged consequent errorsin fluence measure- 
ment to be negligible. 


, 
7 


D. Gamma-Ray 
Counting of ‘e’Au 


We counted the 412-keV gamma rays from the ac- 
tivated gold in the samples with a Ge(Li) detector. 
To eliminate small errors in ampule shapes and 
placements, each sample was rotated at 0.1 rev/s 
while it was counted. 
The 
sample holder was 
mounted 
rigidly 
and 
reproducibly 
by 
a clamp 
around the case of the Ge(Li) detector, as shown in 
Fig. 8. We calibrated the gamma-ray detector by 
placing 
standard 
sources 
in 
the 
same 
sample 
mounting. One standard source* was 152Eu,which 
has, by coincidence, a 412-keV line. A second stan- 
dard source was a few milligrams of activated gold 
wire that LASL radiochemists (Group CNC-11) 
counted in a NaI(Tl) well detector and also with a 


GPurchased from Laboratoire de M6trologie des Rayonnementa 
Ionisants, Commissariats a l’Energie Atomique, Saclay, France. 


dN2 


~ 
“ ‘1NAuU2$ - 
A2N2 
, 
(2) 


where K1 is the ratio of neutron fluence averaged 
over the length of gold wire to that averaged overt he 
volume of lithium hydride (see Appendix B for this 
calculation); a is the thermal capture cross section 
of gold; and AXis the decay constant of lqlAu. The 
solutions 
of 
these equations for 
an irradiation 
starting at t = Oand ending at t = t. are 


Nl(te) 
= NLi-6CJ1[l 
- 
exp(-~lte) 
]l$l/X1 
(3) 


N2(te) 
= K1NAUIJ2[1- 
exp(-k2te) 
]1$/A2 
, 
(4) 


and these give 


ual 
- 
‘Li-6 
1 2 
exp (-A1te) 
Nl (t=) 
= N2(te) 
NAUU2X1 1 - 
exp (-A2te) 
“ 
(5) 


Ge(Li) 
detector. 
Both 
of 
these 
counters 
were 
calibrated with a NBS mixed radionuclide emission 
rate standard. 


E. Calculation of Expected Tritium Activity 


To determine the relation between the predicted 
tritium decay rate dNJdt and the number of ‘98Au 
gamma rays counted, we note that the net produc- 
tion rate of tritium is 


where al is the thermal cross section for ‘Li(n,a), 
AI 
is the decay constant for tritium, and I#Jis the ther- 
mal neutron fluence averaged over the volume of the 
ampule. Likewise, the net production rate for lUIAU 
is 


Fig. 8. 


Rotating sample holder mounted on the Ge(Li) 
detector 
that was used to detect 
the 412-keV 
gamma 
rays from 
test samples 
activated 
by 
thermal 
neutrons. 


We determined N,(tJ, the number of ISSAUatoms 
at the end of the irradiation, by “counting the 412- 
keV gamma rays with a Ge(Li) detector of efficiency 
e. The number of ‘g8Audisintegrations per second in 
the interval At starting at elapsed time t, after the 
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end of the reactor irradiation was determined from 
D, the number of 412-keV gamma rays detected, us- 
ing the relation 


where 
(1) 


(2) 


(3) 


(6) 


T is the measured transmissionof the 412-keV 
gamma rays through the LiH sample and 
quartz wall; 
the 
constant 
/3 is a 
combination 
of 
the 
branching ratio and internal conversion coef- 
ficient in the decay of “Au; 
there are @ 412- 
keV gamma rays per disintegration; 
K, is a coefficient that relates to the size of the 
gold wire used to measure neutron fluence 
(some neutrons entering the gold wire and 
gamma rays leaving the wire are absorbed). 
In Appendix C we calculate KZ,the ratio of the frac- 
tion of gamma rays escaping the gold to the limiting 
value for an infinitely thin wire. Now 


E(tl, At) - 
N2(te) 
exp(-A2t1) 
[1 - 
exp(-X2At) 
] 
, 


(7) 


where we have accounted for ‘BEAudecay during 
elapsed time ti starting at time t. and also for decay 
during the counting interval At. From Eqs. (6) and 
(7), 


D exp(lztl) 
‘z(te~ 
= TCBK2 
[I - exp(-~2At)1 


and substituting for NZ [Eq. 


9 
(8) 


(5)] we get for the 
number of tritium atoms formed in terms of the 
number of 412-keV gamma rays detected, 


N 
= ‘Li-6”lA2’1 
- 
‘m ‘-Alte) 
] 
1 
NAuu2~l [1 - 
exp (-A2te) 
1 


D e.xp(A2t1) 


“ Te8KlK2[1 - 
exp(-A2At) 
J 
“ 
(9) 


The mass of lithium per gram of LiH, g, was deter- 
mined by Union Carbide’s analysis of the material. 
T%e numbers Nu., and NAUof lithium and gold 
atoms are given by 


and 


A 
.—~ 
‘Au 
‘Au 
2’ 


(lo) 


(11) 


where AL1and AA. are the atomic weights (for the 
isotopic mixtures as used) of the lithium and gold, A 
is Avogadro’s number, and f is the isotopic abun- 
dance of ‘Li. The number of tritons formed per unit 
mass of LiH is then 


‘1 
a1A2g 
_=%. 
, 
fD exp(A2t1) 


5 
%1 
“2A1”2 
Te6K1K2 


1- 
exp (-Alte) 


“[l- 
(12) 
e.xp(-12te)1[1 - exp(-A2@ 1 
“ 


To 
compare 
calculated 
and 
measured 
tritium 
counting rates, we determined the decay rate per 
unit mass of LIH: 


AAuu1k2gf [1 - exp(-A1te)1 D em (a2tl) 


~iu2rn2[l - em (-k2te) 1[1 - exp(-a2At)]TeBK1K2 “ 


(13) 


The constants appearing in Eqs. (12) and (13) have 
the values given below. 


Conetant 
value 


196.W 


7.01606 
936*4b 
98.8 + 0.3 b 
1.7781x 10-’ S-1 
2.976x 10-’ s-’ 
0.67075 
0.9447 (large ampulea) 
0.9710 (small ampules) 
3.9611 x 10-’ 
0.9547 
9.42 X 10-’ 


1.005 (small ampules) 
1.013 (large ampules) 
O.aaao 


Reference 


CRC Handbook of Chemist~ 
and Physicu, 56th ed. 
(calculated) 
8 
9 
10 
11 
Union Carbide 
(meaaurad) 
(measured) 
(measured) 
12 
Union Carbide and LASL 
analyses 
Appendix B 
Appendix B 
Appendix C 
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Table IV, which gives the predicted [from Eq. 
(13)] and the measured tritium activities for the 
lithium 
hydride 
samples irradiated by 
thermal 
neutrons, shows that predicted-to-calculated ratios 
have values characteristic of the size of the sample. 
Table V summarizesthis dependence and shows also 
the ratios of predicted and measured volumes of 
evolved 
gas. The 
fact 
that these have similar 
systematic variations suggesta that most of the 
systematic error in Teledyne’s tritium measurement 
is characteristic of the process of gas evolution and 
conversion to water. Using the thermal fluence 
tritium-production measurements to establish stan- 
dards, we normalized the tritium determination 
results from Teledyne for the small and large sam- 
ples by the divisors (from Table V) 1.103 and 1.058, 
respectively. If we note that Teledyne’s 
quoted 


systematic error for the entire proceiw was 5Y0,we 
find these normalizations quite plausible. 
We adduce an additional piece of evidence that 
supports the assumption that the systematic errorin 
the Teledyne measurement is characteristic of the 
sample processing. Teledyne measured the tritium 
decay in sample 305, which contained 0.1594 g of 
LiH, with 0.1 at.% CLi.This was one of the test sam- 
ples irradiated in thermal neutron flux. This sample 
was converted to gas and additional hydrogen carrier 
gas was added in a manner typical of all small am- 
pules. The total volume was then converted to 
water; an aliquot of this sample was converted back 
to hydrogen and counted in a gas proportional 
counter 
with 
a 
resulting 
activity 
of 
347.15dis OS-’o (g H,O)-’ as of 24 June 1976,Another 
aliquot of the water from this sample was certified 


TABLE 
IV 


TRITIUM 
ACTIVITY 
INDUCED 
BY THERMAL 
NEUTRON 
IRIMDIATION 


Sample No. 


Mass ZiH, 


:) 


302 
305 
307 
309 
312 
313 


315 
319 
321 
326 
32a 


0.1384 
0.1594 
0.1655 
0.1598 
0.1496 
0.1682 


0.8789 
0.9269 
1.0484 
0.9291 
0.9232 


l?ritium Counting Rate,” 
nr,l odNr/dt 
[10’ counts ”s-’”(g 
LiH)-’] 


Predicted 
Measured 


2.205 
2.362 
2.165 
2.528 
2.135 
2.362 
2.067 
2.212 
2.271 
2.528 
2.150 
2.345 


2.223 
2.428 


2.180 
2.262 
2.319 
2.462 
2.117 
2.191 
2.283 
2.428 


Ratio 
Measured/ 
Predicted 


1.071 
1.168 
1.106 
1.070 
1.113 
1.091 


av = 1.103 
* 0.015 


1.092 
1.038 
1.062 
1.035 
1.064 


av = 1.058 
+ O.O1O 


‘Ae meaeuredby Teledyne and as predicted from the gold activation. Fractional ‘Li abundance was 
9.42 x 10-’. Reference time was 1200 MST, 26 February 1976. 


,. 


.. 
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TABLE 
V 


RATIOS 
OF MEASURED-TO-PREDICTED 
VALUES 
FOR 
VOLUMES 
AND 
ACTIVITIES 
OF EVOLVED 
GASES 


Average 
Ratio 
Measured-to- 
Sample 
of Volume 
Predicted Tritium 


size 
H2” 
Activityb 


Small 
1.0257 + 0.005 
1.103 * 0.015 
Large 
0.9519 & 0.015 
1.058 + 0,010 


‘Predicted values were obtained from LiH mass and the Union 
Carbide reported value for the number of grame of H, per gram of 
LiH. 


b%edicted valuea are derived from thermal activations (see Table 
Iv). 


by NBS to have an activity of 315.5 + 6 dis. s-’. 
(g H,O)” 
as of 24 June 1976. (The uncertainty of 6 
s-is g-l is the linear sum of 1.3, the random error at 
99% confidence level, and 4.7, the linear sum of con- 
ceivable systematic errore.) The ratio of these two 
measurements is 1.1003, consistent with the correc- 
tion factor (see Table IV) determined by the gold- 
monitored thermal activations of the small samples. 
We have already seen (Table II) that Teledyne and 
NBS determinations of tritium activities in two frac- 
tions of a sample of tritiated water agreewithin 1.5Y. 
(which is expected because both determinations are 
referred to as the same standard sample) so the 
systematic 
error is probably 
characteristic 
of 
Teledyne’s gas handling system. 


IV. 
TRITIUM-PRODUCTION 
RESULTS 


The neutron production from the T(d,n) reaction 
in 
the 
target 
inside 
the 
aLiD 
assembly 
wa5 


monitored by counting the alpha particles from the 
target in a small solid angle at 135° with respect to 
the deuteron beam direction. The total number of 
source neutrons produced was (9.42 + 0.28) x 1016. 
Additional data on the total number of neutrons 
and the change in their energy spectrum as they 
were transported through the ‘LiD assembly were 
acquired by the LASL Radiochemistry Group, CNC- 


11, through the use of various activation foils. These 
resultda appear in Part II; those obtained at the in- 
ner surface provide 
a measure of 
the 
14-MeV 
neutron flux and are compatible with the alpha- 
particle monitor data. 
The results of Teledyne’s analyses of the samples 
of 
lithium 
hydride 
are given in Table 
VI. 
In 
tabulating the Teledyne counting data, we sub- 
tracted from their reported number of counts per 
gram of lithium hydride the same quantity for unex- 
posed samples. The backgrounds measured in ir- 
radiated 
ampules 
containing 
only 
helium 
were 
negligibly small. We applied an additional correc- 
tion because Teledyne used 12.26 yr for the half-life 
of tritium in determining the activity as of 26 
Feburary 1976 from an NBS 
standard dated 3 
September 1961. The apparent best value for the 
tritium half-life is 12.35 yr (Ref. 10), implying that 
Teledyne’s counting data should be multiplied by 
1.006. 
The results for both lithium isotopes are plotted in 
Fig. 9, along with the results of measurements by 
Wyman.1 Although his measurements were made in 
an assembly of the deuteride of natural lithium, his 
results are nearly the same as those in this experi- 
ment. 
The fractional standard deviation for data listed 
in Table VI is 6Y0. 


To evaluate the systematic errorsin these tritium- 
production 
measurements, we assumed that the 
standard deviations in the normalizations (those 
found for the measured-calculated ratios in Table 
V—1.5 and 
1.0% for small and large ampules, 
respectively) 
are also systematic errors in these 
determinations. The standard used to determine the 
counter efficiency for the 412-keV gamma rays for 
gold had an error quoted by the Laboratoire de 
Metrologiesdes Rayonnements Ionisants as +2% at 
%ly. confidence 
level. 
We 
checked 
this source 
against a separate standard, used by radiochemists 
at LASL, whose uncertainty was also quoted as +2Y0 


at the 99% level. 
The 
two standards 
agreed 
to within 


2Y0, and 
we used 
an average 
of the 
two, 
assigning 
a 


2~0 systematic 
error. 
The 
standard 
deviation 
of ayo 


in the neutron-source strength also appears as a 
systematic error. We used the linear sums of the er- 
rors to obtain total systematic errors of 6.5 and 6°% 


for small 
and 
large 
samples, 
respectively. 
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TABLE 
VI 


RESULTS 
OF TRITIUM-PRODUCTION 
MEASUREMENTS 


TritiumDecay Rate’ 
[s-’*(g 
LiE)-’] 
Lacation 
(seeFig. 4) 
Mnms 


Normalised 
f(r)’ 


Teledvne 
and 
110-” atoms.mm* 
Teledyne 
(U26176) 


33.0 
36.67 
46.0 


226.3 
240.0 
240.0 


., 


., 


Li 
Sample 
Ioo@e 
No. 
— 
— 


6 
136 
137 
146 


138 
139 
140 


Lilr 
(6) 
(LA) (:;) 


“Ne@- 
(2/26/76) 


(krrrectad” 
- 
“(Liatom)-’ 


(3/27/76) 
o(source neutron)-’] 


1.107 
1.0544 
1.0429 


29%5 
299.5 
299.5 


-: 
135 


-15 
-3s 
-145 


-30 
-135 
30 


-55 
40 


-4 
64 


–20 
120 


–126 
35 


-30 
30 
136 


–40 
55 


-60 
46 


31.23 
24.66 
34,71 
27.41 
42,69 
33.63 


216.1 
77.24 
0.9464 
1.0213 
1.0584 


0.9963 
o.a631 
1.0939 


0.1409 
0.1218 


0.1432 
0.1367 


0.9054 
0.9128 


0.s432 
0.9343 


201.6 
201.6 
201.5 
227.2 
227.2 
81.21 
81.21 


689.6 
662.1 
832.5 


81.50 
79.01 
87.70 


141 
142 
143 


104 
111 


127.5 
127.6 
125.5 


601.7 
583.3 
868.3 


1067 
1086 
77.64 
76.86 
969 
988 
51.42 
61.13 


106 
110 


7 
163 
1s6 


147 
173 


1717 
1883 


17.17 
17.0 


52.67 
69.33 


1559 
1526 
36.89 
32.97 
61.2 
49.5 


299.6 
289.5 
13.39 
13.22 
12.67 
12.51 
11.36 
11.21 


201.6 
199.5 
47.89 
54.65 
45.33 
61.83 
18.39 
20.63 


149 
166.33 
171 
166.7 
172 
166.7 


101 
570.0 
114 
666.7 


103 
1353 
108 
1342 


‘Isotopic abundances in the test samples were: ‘Li, 95.6034%; T& 99.9058%. MolecuIar weights used were ‘LiH, 7.066; 
7LIH, 8.023. 


bPositive angles are below the parting plane in Fig. 4; negative angles are ahove. 


‘Counting rate of evolved tritium. 


%itium 
activities in %ii 
samples were corrected for a background measured in unirradiated samples of 3.78 + 0.52 s-’ 
. (g LiH)-i 
(see Table III); ‘LiH samples had negligible background. 


Teledyne 
counting data were corrected because they were based on a tritium half-life of 12.26 yr instead of the presently 
accepted 12.35 yr; this affecta the NBS standard used (dated 9/3/61). 


The specific tritium production, f(r), (equals 4+ 
times number of tritons formed per atom of lithium for one source 
neutron). 


0.7892 
0.9068 
0.8788 


0.1465 
0.1211 


0.1340 
0.1303 


184.6 
162.9 
162.9 


166.8 
154.1 
164.1 


25.31 
25.03 
25.03 


127.5 
125.5 
125.5 


77.68 
75.66 


668.2 
562.9 
514.0 
611.0 


30.96 
26.23 


49.5 
61.6 
1349 
1338 
1224 
1215 
29.97 
32.20 


. 
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Fig. 9. 
ZYitium productwn 
in a sphere of ‘LiD plotted 
as f(r) = 41r? times tlw number of tritow 
per 
atom of lithium for one 14-MeV neutron at the 
center 
of 
the 
assembly. 
The 
triangles 
are 
measurements 
by Wymard of tritium 
produc- 
tion in 7Li metal embedded 
in a sphere of LiD. 
The circles are our ‘Li production 
data; 
the 
squares, 
‘Li data. 
The curves 
are drawn 
to 
guide the eye. 
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APPENDIX 
A 


TRITIUM 
DECAY 
RATES 
IN SAMPLES 
OF UNIRRADIATED 
7LiD 


Tritium decay rates in samples of 7LiD prepared 
TABLE A-I 
by Union Carbide Corporation are listed in Table 
A-I. 
BACKGROUNDS 
IN ‘IAD SAMPLESa 


Sample 
No. 


65 
69 
107 
108 
109 
110 
111 


Mass LiD 


(!?0 


0.1733 
0.20% 
1.1078 
1.1397 
1.0285 
0.9681 
1.0162 


Activity 
[10’ counts ”s-’” (g LID)-’] 


1.638 
“ 1.138 
1.317 
1.317 
1.299 
1.242 
1.221 


av 
= 1.310 + 0.060 


“Aeof 1 July 1975. The atomic abundance of %i was 1.047 x 10-4. 
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CALCULATION 


To find a value for the thermal 


APPENDIX 
B 


OF VOLUME-AVERAGED 


neutron fluence 
averaged over the volume of lithium hydride in a 
quartz ampule, we must take account of the at- 
tenuation of flux as it passes through the lithium 
hydride itself. These types of calculations have been 
done for purely absorbing media by Case, de Hoff- 
mann, and Placzeki’ for several different geometries. 
For a sphere of radius 
a 
immersed in an isotropic 
bath of monoenergetic neutrons, the average flux 
throughout the volume relative to that far away is 
given by 


PO=+; 
(?‘< ‘) 
‘ 


where 1?is the absorption mean free path. For a 
sphere of lithium hydride [(9.42x 10-2)% CLi]with a 
= 8.0 mm and.4 = 257.9 mm, a/,C= 0.0310 and PO= 
0.9767. Thus the fluence averaged over the sphere is 
only a few per cent different from the fluence far 
away. This difference is an upper limit of the correc- 
tion necessary because the gold was distributed 
along a diameter and not throughout the volume. 
If 
o(r) 
represents the 
radial 
distribution 
of 
neutrons in the sphere, then the average fluence in- 
teracting with the lithium hydride is given by 


a 


TV= 31a3J 


f$(r)r2 
dr 
, 


0 


NEUTRON 
FLUENCE 


and the fluence averaged along a radius is given by 


a 


7=- 
Ila J_ 


$(r) 
dr 
. 


-o 


With the gold wire positioned along a diameter of 
the ampule, we determined & 
by counting the 
number of 412-keV gamma rays. 
We used the code DTF-IV (Ref. 15) to estimate 
O,(r) in one-group S. calculations with n = 8, aasum- 
ing isotropic scattering in the center of mass. 
Calculations were performed using both purely ab- 
sorbing media and with several different estimates 
of the scattering cross section for lithium hydride 
molecules. Neutron interactions in the quartz am- 
pules were also included in the calculations. We 
calculate that the ratio of the averaged fluxes, 


is 1.005 + 0.002 for the small ampules and 1.013 + 
0.005 for the large ones. This relation was used to 
compute the amount of tritium produced in the two 
different size ampules. 


APPENDIX 
C 


NEUTRON 
AND 
GAMMA-RAY 
ABSORPTION 
IN 
GOLD 
WIRE 


.* 


.. 


. 


0 


A 0.05-mm-diam gold wire in the form of a helix 
in a fluence lower than that at the outside of the 
was positioned along the diameter of the ampule to 
wire. Furthermore, some of the 412-keV gamma rays 
measure the thermal neutron fluence within the test 
originating in the interior of the wire fail to escape. 
samples of LiH. Gold has a sufficiently large neutron 
These effects are smallerthan 5%, and we calculated 
capture cross section that the interior of the wire is 
the correction for them. 
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.. 


We assume that the pitch of the gold helix is great 
enough—perhaps 50 wire diameters—that we can 
make the calculation for an infinite straight wire. 
The calculations for neutrons entering and gamma 
rays leaving are identical. We referto the calculation 
of average fluence by K. M. Case, F. de Hoffmann, 
G. Placzek14for a purely absorbing medium. Their 
calculation is appropriate for gold, which has a cross 
section for absorption 10 times greater than that for 
scattering. For radii small compared to the mean 
free path, they find for the escape or entrance 
probability, 


——.——— 
———— 
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PART 
II 


CALCULATIONS 
OF TRITIUM 
PRODUCTION 
AND 
RADIOCHEMICAL 
ACTIVATION 
IN A SPHERE 
OF 
aLiD IRRADIATED 
BY 
14-MeV NEUTRONS 


by 


Jon M. Wallace 


. 


ABSTRACT 


The SPCCWICproduction of tritium in ‘LiH and 7LiH embedded ampules and 


the activation of radlochemical 
detector foils of %3c, 89Y,‘Zr, 
10%, 
19%.879 
10gIr.,2,,lnAu, ‘a’U, and 28SUplaced at various positions in a 600-mm-diam 
‘LiD 
sphere irradiated 
by a central source of 14-MeV neutrons have been 
calculated 
and 
compared 
with 
experimental 
data. 
One- 
and 
three- 
dimensional Monte Carlo and S. neutron-transport 
calculations 
were per- 
formed. 
The 
most reliable 
(three-dimensional 
Monte 
Carlo) 
calculation 
generally 
agrees with both the tritium-production 
and the radiochemical- 
activation data. The existing discrepancies between calculation and experi- 
ment 
appear 
largely 
attributable 
to 
uncertainties 
in 
production and radiochemical-activation 
cross sections. 
some 
tritium- 


1. INTRODUCTION 


We 
present 
the 
results of 
calculations in a 600-mm-diam 
neutron-transport 
sphere of OLiDir- 
radiated 
by 
a central 
14-MeV 
neutron source. 
Tritium production from lithium in embedded aLiH 
and 7L1Hampules and the neutron-induced activa- 
tion of radiochemical detector foils placed at various 
positions in the sphere were obtained and compared 
with 
experimental data. 
The 
‘LiD 
sphere, the 
neutron 
source, 
and 
the 
tritium-production 
measurements have been described in detail in Part 
I of this report. Additional diagnostic information 
was obtained from six radiochemical foil packets 
placed 
at the interfaces of the concentric ‘LiD 
spherical assembly. The foil materialswere‘“SC, ‘SY, 
‘Zr, 
l’~m, 
“’Ir.,,, ‘%.,,,, 
“’Au, “W, and ‘“”U. The 
foils were analyzed by radiochemical methods after 
the irradiation. 
The comparison of our calculations with the in- 
tegral experiment provides a simultaneous test of 


the neutron-transport codes, 
the neutron cross- 
section data relevant to transport through ‘LiD, the 
tritium-production 
cross-section 
data, 
and 
the 
radiochemical-activation cross-section data. 


II. 
CALCULATIONS 


Four different calculations of this experiment will 
be discussed. They vary in reliability, and their in- 
tercomparison is instructive because it points out 
the shortcomings of the less accurate calculations. 
Each 
calculation 
consisted of two parts: 
(1) 
a 
neutron-transport computation to produce fluences, 
that is time-integrated fluxes, at various detector 
positions, 
and 
(2) 
an 
evaluation 
of 
the 
radiochemical-activation integrals and, for Calcula- 
tions I-III, the tritium-production integrals for com- 
1 


parison with the experimental values. In Calculation 
IV, the tritium production was obtained directly in 
the neutron-transport computation. The ‘Li and 7Li 
I 
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I 


cross-section data used in both the transport and 
tritium-production computations were ENDF/B-IIL 
The 1967 United KingdomfLASL evaluation of the 
deuterium cross sections was employed in the tran- 
sport computation. 
The radiochemical-activation 
cross sections used werefrom the LASL TD-Division 
Dosimetry Library.* The (n,f) and (n,~) cross sec- 
tions in this library werenormalized using BIG 10** 
data. There is a 3% normalization uncertainty in our 
results that is due to the uncertainty in the number 
of 14-MeV source neutrons, as described in Part I of 
this report. The calculations will now be discussed in 
order of increasing reliability. 


A. 
Calculation 
I—lD, 
1l-Group, 
S,,P1 


In this calculation the neutron fluences at the re- 
quired radial distances from the source were com- 
puted using a one-dimensional neutron transport 
code with n-group, 
St,Pl neutron transport. The 
group structure is given in Table I. The source 
neutrons were distributed between 13.5 and 17.0 
MeV according to an appropriate weighting func- 
tion. 
The 
radiochemical-activation 
and tritium- 
production integrals were computed with cross sec- 
tions collapsed into the n-group structure using the 
same weighting function. 


B. 
Calculation 
II—ID, 
21-Group, S8,P1 


This calculation is similar to Calculation I, except 
smaller energy and angle bins are used in a 21-group, 
S,P, transport scheme. The energy group structure is 
given in Table II. In this computation, group 1 is 
empty. 
The 
source 
neutrons 
were distributed 
between 
13.5 and 
15.0 MeV 
according 
to 
the 
weighting function used in Calculation L As in 


*The Sc(n,-y) cross sections are based on ENDF/B-IV, Dosimetry 
Library. The Y(n,y) cross sections were taken from an evaluation 
by E. D. Arthur in “Applied Nuclear Data Research and Develop- 
merit, ”LASL report LA-6971-PR (September 1977) (compiled by 
C. I. Baxman and P. G. Young), p. 3. The Tm(n,y) cross sections 
are from a private file of J. S. Hendricks. The (n,2n) cross sec- 
tions, except for ‘“U, were taken from evaluations based on data 
in B. P. Bayhurat, Phys. Rev. C 12, 451 (1975). All other activa- 
tion cross sections were based on ENDL (Llvermore Evaluated 
Nuclear Data) 4/76. 


**BIG 10 is a critical assembly at LASL. 


TABLE 
I 


ENERGY 
GROUP 
STRUCTURE 
FOR 
CALCULATION 
I 


Group 
Number 


JL.ow. 
(MeV) 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 


13.5 
10.0 
7.79 
3.68 
2.232 
0.500 
0.184 
0.0248 
0.00335 
0.000167 
0.139 x 10-9 
EM-X= 17.0 


Calculation I, this weighting function was also used 
to 
collapse 
the 
radiochemical-activation 
and 
tritium-production cross sections into the 21-group 
structure. 


C. Calculation 
111—lD, Monte Carlo 


In this calculation the neutron fluences at the re- 
quired radial distances from the source were com- 
puted using the three-dimensional Monte 
Carlo 
neutron-transport code MCN.aIsA one-dimensional 
mockup of the spherical assembly, similar to that in 
Calculations I and II, was used. The fluences, which 
were surface integrated, were binned into a 75-group 
structure (Table 
III). The point-source neutrons 
were distributed isotropically with energy 14.1MeV. 
Radiochemical-activation 
and tritium-production 
integrals were calculated using precollapsed cross 
sections and the 75-group histogram fluences. The 
fluences in the high-energy range relevant to the 
(n,2n) activations were increased by 2% to simulate 
the effect of the observed source anisotropy in the 
direction of the radiochemical detectors. This is our 
most reliable one-dimensional calculation. The rele- 
vant geometric parameters for Calculations I-III are 
given in Table IV. 


., 


., 


. 


, 
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TABLE 
II 


ENERGY 
GROUP 
STRUCTURE 


FOR 
CALCULATION 
II 


Group 
Number 


&ower 
(MeV) 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 


15.0 
13.6 
12.0 
10.0 
7.79 
6.07 
3.68 
2.865 
2.232 
1.783 
1.363 
0.823 
0.600 
0.303 
0.184 
0.0676 
0.0248 
0.00912 
0.00336 
0.000167 
0.139x 10-’ 
EM.. = 17.0 


TABLE 
111 


ENERGY 
GROUP 
BOUNDARIES 


FOR 
BINNING 
IN 
CALCULATIONS 
111AND 
IV 


16.22MeV’ 
15.12” 
15.00” 
14.87” 
14.75’ 
14.62” 
14.50” 
14.37” 
14.25’ 
14.llb 
14.12’ 
14.05b 
14.00 
13.87 
13.75 
13.62 
13.50 
13.37 
13.25 
13.12 
13.00 
12.83 
12.67 
12.50 


12.33 
12.16 
12.00 
11.75 
11.50 
11.25 
11.00 
10.75 
10.60 
10.25 
10.00 
9.75 
9.50 
9.25 
9.00 
8.75 
8.50 
8.16 
7.79 
7.40 
7.00 
6.75 


6.60 
6.28 
6.07 
5.78 
5.50 
5.26 
5.00 
4.75 
4.50 
4.09 
3.68 
3.27 
2.865 
2.55 
2.232 
2.000 
1.738 
1.550 
1.353 
1.090 
0.823 
0.660 


0.500 
0.400 
0.303 
0.244 
0.184 
0.125 
0.0676 
0.0460 
0.0248 
0.0170 
0.00912 
0.00600 
0.00335 
0.001235 
0.000454 
0.000167 
o.m14 
O.0000 


‘Energy boundaries for Calculation lV, 


Wnergy boundaries for Calculation III, 


TABLE 
IV 


GEOMETRICAL 
DATA 
FOR 
CALCULATIONS 
I-III 


Outside 
Radius 
Density 


(mm) 
Material 
(g/cm’) 


22.2 
50.0 
51.0 
76.15 
77.15 
126.0 
127.0 
200.0 
201.0 
300.0 


0.0 
0.753 
0.0 
0.763 
0.0 
0.753 
0.0 
0.751 
0.0 
0.755 
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D. 
Calculation 
IV—3D, 
Monte Carlo 


The neutron fluences in this calculation also were 
obtained using MCN. The calculation incorporated 
a 
simplified 
three-dimensional 
mockup 
of 
the 
spherical assembly, including small spherical ‘LiH 
and 7LiH regions to simulate the tritium collection 
ampules. The geometric data are given in Table V. 
The 
geometry is illustrated in Fig. 
1. Tritium 
production was calculated on-line by the neutron- 
transport 
code. 
Neutron 
fluences 
for 
the 
radiochemical-activation 
integrals were tallied at 
the points where the detectors were located for 
neutron energies >7.40 MeV. 
Surface integrated 
fluences, as in the one-dimensional calculations, 
were used for energies E z 7.40 MeV. This hybrid 
treatment was employed because the MCN point- 
fluence tally is very time consuming. The tally 
should not be necessary at lower energies because 
the multiple-scattering process, which produces the 


TABLE 
V 


GEOMETRICAL 
DATA 
FOR 
CALCULATION 
IV’ 


Density 
Material 
(g/cms) 


aLi.s~6~7Li.0~~~ 
D 
0.7425 
Stainless steel (SS) 
9.03 
LiH 


Distance of Center 
From Source 
(mm) 


0,300 


Radius 
of LiH Spheres 


Experiment 
(mm) 
Calculation 
(mm) 


50.5 
76.65 


126.5 
200.5 


300.0 


5.0 
5.0 


9.0 
9.0 


9.0 


5.0 
6.0 ‘LiH 
5.5 lLiH 
9.0 
9.0 ‘LiH 
12.07LiH 
15.3 


“Source volume: spherical region of radius 2.0 mm. 
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Fig. 1. 
The 
geornetqy 
of the 
‘LiD 
sphere 
used 
in 
Calculation 
IV. 
The 
small 
embedded 
LiH 
spheres are not shown, 
but are located at the 
positions 
illustrated 
in Fig. 4 of Part I of this 
report. 
The radii 
of the LiH spheres 
in the 
calculation 
are in some cases larger than the 
radii of the ampules 
used in the ezpen”ment 
(Table 
V). 


lower energy component of the neutron spectra, 
tends to produce an isotropic distribution as the 
neutrons are downscattered. For this calculation, 
the fluences were binned into an 83-group structure 
(Table III), which consisted essentially of the group 
structure of Calculation III with eight additional 
groups at the upper energy end. The calculation 
employed 
an anisotropic energy-angle correlated 
source occupying a spherical volume, which closely 
simulated the experimental source (see Fig. 3 of Part 
I of this report). The energy-angle correlation is that 
produced by a 240-keV deuteron beam on a tritium 
target.* The 
radiochemical activations were ob- 
tained just as in Calculation III, using the histogram 
fluences. Of the four calculations, this one is the 


*For this calculation, this correlation is identical to that for the 
300-keV beam actually used in the experiment. 


most precise simulation of the experiment. Com- 
parison of the results of this calculation with the ex- 
perimental data provides the best test of how well we 
understand the processes under consideration. 


III. 
TRITIUM 
PRODUCTION 


The experimental tritium-production results at 
each radius have been averaged for comparison with 
the calculations. This average is the proper com- 
parison for Calculations I-III, which are spherically 
symmetric, and provides better statistics for com- 
parison 
with 
Calculation 
IV. 
In 
addition, 
the 
tritium-production results are compared ampule by 
ampule with the results of Calculation IV. 
All comparisons are made in terms of the ratio of 
the experimental results to the calculated results. 
The quoted uncertainties for this ratio are from the 
experimental uncertainties only. 


A. Tritium Production From ‘Li 


The 
material contained 
in the ampules was 
‘Li.SOw,8‘Li.,m,2 H. So the quantity observed was 


‘~i-’ 
= 4nr2[00’’’8~0,(E)uL,L7(E)E) 
‘E 


+ 0.000942 
J 
1 


I$i(E)uLi_6(E) 
dE 
s 


where (hi.1 is the reaction cross section for 7Li(n,n’a) 
and ~IJ.@is the reaction cross section for ‘%i(n,a). 


Here O, is the neutron fluence spatially integrated 
over the ampule volume. The amount of tritium 
produced from ‘Li is never more than 0.3% of the 
total. The ‘L1 tritium-production reaction threshold 
is at 2.8 MeV; thus these data probe neutron tran- 
sport only in the energy range above this threshold. 
The results of the calculations are presented in 
Tables VI and VII and shown in Fig. 2. The various 
calculations differ from each other by as much as 


‘zs~. 
at the 300.O-mmradius (at the outside of the 
LID sphere) with typically 10% differences at the 
other positions. All calculations give overpredictions 
at the inner three positions. Calculation IV gives the 
best overall agreement with experiment, the dis- 
crepancies 
being 
generally smaller at 
positions 
further away from the source. These results are 
similar to those of a recent analysis of the Wyman 
experiment.4-6In Sees. lV.A and IV.B we discuss 
(n,2n) and (n,f) activations of ‘W, 
in which there is 
generally good agreement between calculation and 
experiment. The cross sections for these two reac- 
tions are qualitatively similar to the 7Li tritium- 
production 
cross section and are thought to be 
relatively well known. Hence it is not unreasonable 
to attribute the 7Li discrepancies to the tritium- 
production cross section used. A 15% reduction of 
the cross section in the 14-MeV region would bring 


TABLE 
VI 


RATIO 
OF OBSERVED-TO-CALCULATED 
TRITIUM 
PRODUCI’ION 
FROM 
‘Li 


Distance From 
Soum% 
Average Experimental Result 
(mm) 
(10-” atoms ”mm’/Li atom) 


50.5 
31.09 
76.65 
30.10 
126.5 
25.12 
200.5 
19.46 
300.0 
11.29 


‘Ratio of experiment to calculation. 


Calculation’ 


I 
II 
III 
IV 
—. 
—. 


0.783 
0.836 
0.873 
0.883 
0.811 
0.838 
0.882 
0.861 
0.834 
0.809 
0.860 
0.890 
0.989 
0.883 
0.926 
0.918 
1.266 
1.024 
1.015 
0.986 


TABLE 
VII 


CALCULATION 
IV AND 
OBSERVED 
TRITIUM 
PRODUCTION 
FROM 
‘Li 


Sample 
No. 


109 
103 
114 
101 
171 
149 
172 
173 
147 
183 


Distance 
From Source 
(mm) 
Observation 
(10-” atoms”mm’/Li atom) 
Calculation IV 
(10-” atoms ”mm’fLi atom) 


50.5 


76.65 


126.5 


200.5 


300.0 


32.20 
29.97 
29.23 
30.96 
25.03 
25.31 
25.03 
20.53 
18.39 
11.36 
11.21“ 


35.6 
34.8 
35.5 
34.4 
223.3 
29.5 
26.9 
21.1 
21.3 
12.0 
10.9 


the calculation into significantly better agreement 
with the experiment. 


B. Tritium Production 
From ‘Li 


The 
material contained 
in the 
ampules was 
7Li.,,,M,‘L1.,,,,,, H; so the observed quantity was 
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TRITIUM 
PROWCTION FRCM 
h 


I 


..,~ 
0 
!s0 
1001s02002!s0300 


RADIUS (mm) 


Fig. 2. 
Ratio 
of 
observed-to-calculated 
tritium 
production 
from ‘Li us radius. 
The error bars 
are from 
random 
experimental 
uncertainties 
only. There is an additional 
normalization 
un- 
certainty 
of -3%. 


‘~-’ 
= 4nr2k-044’’’~$,’E)UL7(E)E) 
“ 


+ 0.955034 
J 
1 
I$i(E)oLi_6 
W 
“ 
. 


The 7Licontributes no more than -5”A of the tritium 
produced in these ampules, with the largest frac- 
tional contributions generally occurring in the am- 
pules closest to the source. The cross section for 
tritium production 
from 
‘Li 
is reasonably well 
known. It has a general decrease with energy, with a 
resonance superimposed at 240-keV neutron energy. 
Most of the tritium production arises from neutron 
energies around and below this resonance energy. 
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TABLE 
VIII 


RATIO 
OF OBSERVED-TO-CALCULATED 
TRITIUM 
PRODUCTION 
FROM 
OLi 


Calculation” 
Distance 
From Source 
Average Experimental Results 
(mm) 
(10-’” at~m~ .mm’/Li a~m) 
I 
II 
III 
Iv 
—— 
—— 


50.5 
34.48 
0.877 
1.020 
1.089 
1.039 
76.65 
51.28 
0.840 
0.953 
1.029 
0.863 
126.5 
82.74 
1.006 
1.086 
1.187 
1.087 
200.5 
79.89 
1.071 
1.074 
1,179 
1.014 
300.0 
28,57 
1.W9 
1.213 
1.441 
0.963 


‘Ratio of experiment to calculation. 


bAmpule density p = 0.3041g/cm’. 


TABLE 
IX 


CALCULATION 
IV AND 
OBSERVED 
TRITIUM 
PRODUCTION 
FROM 
‘Li 


Distance 
Sample 
From Sourca 
Observation 
Calculation IV 
No. 
(mm) 
(10-’4 atoms .mm’/Li atom) 
(lO-U atoms ”mm2/Liatom) 


110 
50.6 
106 
111 
76.65 
104 
143 
126.5 
141 
142 
138 
200.5 
139 
140 
136 
300.0 
137 
145 


32.97 
35.99 
51.13 
51.42 
87.70 
81.50 
79.01 
77.24 
81.21 
81.21 
24.66 
27.41 
33.63 


33.0 
33.4 
59.1 
59.7 
75.5 
79.3 
73.5 
81.5 
76.0 
78.7 
30.3 
30.6 
243.1 


. 


Resonance self-shielding is significant here, because 
in Fig. 3. The various calculations differ from each 
. 
strong contributions come from the region of the 
other by 15 to 25%. As for 7LI,Calculation IV gives 
resulting dip in the neutron spectra. Comparison 
the best overall agreement with experiment. An 
with the ‘Li 
tritium-production data provides a 
anomalously large discrepancy occurs at the radial 
much 
more 
stringent test of 
our calculational 
source distance of 76.65 mm, for reasons unknown. 
capability than comparison with the 7Li data. 
The 
three-dimensional 
calculation 
gives 
a 
The results of the various 
calculations are pre- 
significantly 
different radial dependence 
of the 
sented in Tables VIII and IX and shown graphically 
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Fig. 3. 


Ratio 
of 
observed-to-calculated 
tritium 
production 
from 6Li us radius. 
The error bars 
are from 
random 
experimental 
uncertainties 
only. There is an additiorud 
normalization 
un- 
certainty 
of -370. 


tritium production at this point than do the one- 
dimensional calculations. At all other radii, Calcula- 
tion IV is in rather good agreement with the experi- 
ment. Tritium production from ‘Li was quite sen- 
sitive to small changes in the device geometry, 
probably caused, at least partially, by resonance 
self-shielding. For example, there is a surprisingly 
large sensitivity just to changes in the density of the 
small LiH collection spheres. Changing the density 
from 0.30 to 0.50 glcma induces changes in the 
tritium production by as much as 10% in the larger 
LiH collection spheres. Given this high sensitivity 
together with the simplifications used in the Monte 
Carlo simulation of the experiment, we cannot 
reasonably expect much closer agreement between 
calculation and experiment than was obtained. 


IV. 
RADIOCHEMICAL 
ACTIVATION 


Three classes of activation reactions were con- 
sidered: (n,2n), (n,f), and (n,~). The (n,n’) activa- 
tion of I’iIr was also observed. The targets were “SC, 
““Y, ‘Zr, 
‘a’Tm, ‘s’Ir.,,, ‘%.,,,, 
“’Au, ‘*KU,and “W. 


The foil packets were placed in the plane perpen- 
dicular to the horizontal deuteron beam and 60° up 
from the horizontal axis, @ = 60° and 6 = 90° in the 
coordinate system of Part L The foil distances from 
the source were 22.2, 60.0, 76.15, 126.0, 200.0, and 
300.0 mm. The foils were sufficiently thin that they 
did not significantly effect the neutron flux in the 
‘LiD assembly. 
The observed quantities are 


where Ujis the activation reaction cross section for 
reaction j and 41 is the neutron fluence at the foil 
position i. The data for this part of the experiment’ 
are given in Table X. As in the presentation oft he 
tritium-production results, the results of the activa- 
tion calculations are given as the ratio of the experi- 
ment to calculation and the quoted errorsare due to 
the experimental uncertainties only. 


A. The (n,2n) Activations 


The (n,2n) activations of six different targets were 
observed. The (n,2n) reactions considered here have 
threshold energies ranging from -8 to 12 MeV, ex- 
cept for ‘“aU, which has a threshold of 6.17 MeV. 
These activations consequently result largely from 
the unattenuated 14-MeV component of the neutron 
spectra at the various radii. Comparison of the 
calculations with the measurements provides a test 
of 
how well the transport codes treat 14-MeV 
neutron attenuation. 
In Table 
XI, we give the “Y(n,2n) 
activation 
results for the four calculations, which are typical of 
all the (n,2n) results. Table XII gives the results of 
Calculation IV for all activations. These results are 
shown in Figs. 4-10. From Table XI, we see that 
there is generally good agreement between the ex- 
periment and the two Monte Carlo calculations. The 
S. calculations give less satisfactory agreement with 
the experiment. This is due in part to streaming ef- 
fects in the S“ method, which does not permit 
precisely radial transport. The large discrepancies at 
the smallest radial source distance (22.2 mm), which 
occur in all the calculations, are thought to come 
from 
a slightly off-center 
source. The 
neutron 


., 


., 


. 


. 
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Reaction 


8eY(n,2n)ooY 
‘Zr(n,2n)0sZr 
1’Wm(n,2n)1~m 
W.r(n,2n) l% 
‘nAu(n, 2n)‘WAU 
‘WJ(n,2n)n7U 
Wr(n,2n)lwIr 
+ Wr(n,~)iwIr 
lw~(n,n,)lwrnfi 


2’W(n,f)wMo 
“’U(n,f)wMo 
‘6Sc(n,~)%3c 
‘“Y(n,~)wY 
‘e~m(n,~) ‘7’Tm 
‘WAu(n,~)lwAu 
‘aW(n,~)gaOU 


TABLE 
X 


RADIOCHEMICAL-ACTIVATION 
DATA” 
(10’*) 


Distance From Source 


22.2mm 
50.0 mm 
76.15 
mm 
126.0mm 
200.0 mm 
300.0mm 


1665 


3916 
1540 


2285 


4.592 


72.35 
46.11 


219.2 
145.0 
528.6 
534.7 
578.3 
295.9 
401.1 


61 + 15% 
829 
369 
3.508 
2.029 
59.80 
50.63 
29.75 


76.25 
49.78 
197.4 
197.6 
214.8 
118.7 
166.8 


33 + 25% 
374 
150 
2.492 
1.310 
44.75 
36.44 
20.87 


18.02 
11.44 
55.00 
51.15 
56.88 
37.13 
52.30 


15 + so~o 


142 
46.7 
1.363 
0.8887 
24.38 
19.81 
11.14 


3.999 
2.131+4% 
13.62 
12.70 
13.87 
10.88 
16.20 


8 + 100% 
46.0 
13.0k 4% 
0.5511 
0.2562 
10.04 
8.126 
4.476 


0.6924 
0.3561+ g~o 
2.935 
2.660 
+ 
7yo 


2.994 


2.733 
+ 
5’%. 


2.97 
+ 
7’%. 


8.62+4% 
3.26+ 5% 
0.0971+ ‘7~o 
0.03483 
1.974 * 10% 
1.068 
0.5841 


Wnless otherwise specified, the precision of the data is +3°A. 


TABLE 
XI 


RATIO 
OF OBSERVED-TO-CALCULATED 
“Y(n,2n)wY 
ACTIVATIONS 


Distance From Source 


Calculation 
22.2mm 
50.0mm 
76.15mm 
126.0mm 
200.0mm 
300.0mm 


I 
0.752 
0.740 
0.726 
0.740 
0.803 
0.934 
II 
0.936 
0.926 
0.883 
0.847 
0.848 
0.924 
III 
1.008 
1.067 
0.993 
0.964 
0.938 
1.048 
Iv 
1.026 
1.059 
1.041 
0.983 
0.964 
1.079 


. 


spectrum at this point consists mainly of the unat- 
tenuated 14-MeV component. There is also a low- 
energy component from backscattered neutrons, but 
this contributed negligibly to the (n,2n) activations. 
Moreover, the (n,2n) activation cross sections are 
generally kncpvnto within a few per cent in the re- 
quired energy range. Hence, we disregard the dis- 
crepancies at 22.2 mm. At the other radii, Calcula- 


tion IV seems to be in generally good agreement with 
the experiment. Moat of the remaining discrepancies 
can plausibly be attributed to errors in the activa- 
tion cross sections. However, the results do suggest 
that there may be some small errorsin the total cross 
sections for ‘Li or deuterium in the 14-MeV region. 
This is surmised because the observed-to-calculated 
ratio for a given activation generally decreases with 
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TABLE 
XII 


RATIO 
OF OBSERVED-TO-CALCULATED 
RADIOCHEMICAL 
ACTIVATIONS 
(CALCULATION 
IV) 


Distance From Source 


Reaction 
22.2mm 
50.0mm 
76.15mm 


‘W(n,2n)”Y 
‘Zr(n,2n)”Zr 
‘eWm(n, 2n)’% 
“iIr(n,2n) ‘WIr 
‘oTAu(n,2n)lnAu 
‘a1U(n,2n)a17U 
W.r(n,2n) lo% 
+ W.r(n,~)iwsIr 
lw~(n,n, )Iwmfi 


‘“’U(n,f)”Mo 
‘WJ(n,f)WMo 
“Sc(n,y)’”Sc 
lsY(n,~)mY 
W(’m(n,7)i7”Tm 
‘s7Au(n,~)lWAu 
2uU(n,7)210U 


1,256 
1.035 
0.983 
1.004 
0.999 
1.135 
0.979 
1.136 
1.102 
1.061 


0.744 
1.030 
1.163 
0.992 
1.054 
0.971 
1.015 
0.962 
1.026 
1.059 
1,438 
1.243 


1.1 


m 


1 


II 


I 


I 
f 
09 


0.8 


to.,o~ 


RADIUS (mm) 


Fig. 4. 


1.039 
0.992 
0.997 
0.984 
0.974 
1.057 
1.062 


0.793 
1.001 
0.999 
1.077 
1.028 
0.943 
1.041 
1.189 


E 


a 
9 
$ 


[ 


126.0 
mm 
200.0 mm 
300.0mm 


0.984 
0.941 
1.003 
0.921 
0.939 
1.032 
0.994 


0.908 
1.004 
0.998 
1.084 
1.011 
0.878 
0.983 
1.129 


1.020 
0.849 
0.992 
0.916 
0.925 
1.058 
1.007 


1.463 
0.989 
0.995 
1.085 
0.978 
0.860 
0.964 
1.107 


0.966 
0,819 
0.973 
0.875 
0.919 
1.085 
1.046 


1.047 
1.121 
1.438 
0.991 
1.404 
1.079 
1.166 
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‘zr( 
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Fig. 5. 


Ratio 
of observed-to-calculated 
8’Y(n,2n)8’Y 
Ratio 
of observed-to-~alculated 
‘0Zr(n,2n)”Zr 


., 


. . 


. 


. 


activation 
us radius. 
The error bars are from 
activation 
us radius. 
The error bars are from 
experimental 
uncertainties 
only. 
The 
experimental 
uncertainties 
only. 
The 
threshold 
energy is 11.60 MeV. 
threshold 
energy is 12.12 Me V. 
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Fig. 6. 
Ratio 
of 
observed-to-calculated 
1eeTm(n,2n)1e8Tm activation 
us radius. 
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error bars are from experimental 
uncertainties 
only. The threshold 
energy h 8.11 MeV. 
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Fig. 7. 
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threshold 
energy is 8.17 Me V. 
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Fig. 8. 


Ratio 
of 
observed-to-calculated 
1WAu(n,2n)’WAu activation 
us radius. The error 
bars are from experimental 
uncertainties 
only. 
The threshold 
energy is 8.12 MeV. 
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increasing radial source distance. This behavior is 
particularly pronounced for the highest threshold 
target, ‘OZr(threshold energy 12.12MeV), and is es- 
sentially absent from the “U(n,2n) 
results, which 
have a significantly larger contribution from down- 
scattered neutrons at the larger radii. Optimum 
agreement between the (n,2n) data and calculations 
is obtained by an upward resealing of all ‘Li or 
deuterium cross sections by 4 and 3%, respectively, 
in the 14-MeV region. The production of ‘921ris due 
to 
both 
(n,2n) 
and 
(n, ~) 
reactions. 
The 
leoIr(n,2n)lBiIrreaction contributes more than 70% to 
the total, except at 200 mm where the contribution 
is -56°\0. 


B. The (n,f) Activations 


The 
(n)f) 
activations 
of ‘“U 
and ’80U were 
observed. The mu(n,fl 
reaction has an effective 


threshold of -1 
MeV with most of the activation 
resulting from neutrons above 5 MeV. The reaction 
cross section is thus similar to that of the ‘80U(n,2n) 
and 7Li(n,n’a) reactions. The results of Calculation 
IV are given in Table XII and shown in Figs. 11 and 
12. Disregarding the 22.2-mm data, Calculation IV 
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is in satisfactory agreement with experiment, except 
possibly at the 300-mm radius at the outside of the 
LiD sphere. The reason for this discrepancy is not 
clear, because the fission cross section is thought to 
be known to within a few per cent. 
The 
28VJ(n,f) activation 
is due to 
the entire 
neutron energy spectrum present in this experiment. 
It is the first radiochemical activation discussed so 
far that probes the lower energy region of the 
spectrum. Calculation IV and the experiment are in 
satisfactory agreement everywhere. 
The two Monte Carlo calculations give (n,f) ac- 
tivations that agree everywhere within 1%. The two 
S. calculations are not in such close agreement with 
each other or with the Monte Carlo results. Calcula- 
tion II, the most reliable S“ computation, 
is in 
somewhat better agreement with the experiment 
than Calculation I is. The ‘a6U(n,fl results for the 
various calculations are given in Table X111.The 
2a’U(n,f) results are similar. 


C. The (n,y) 
Activations 


The (n,y) activations of five different target nuclei 
were measured in this experiment. The (n,~) cross 
sections for these targets generally rise rapidly with 
decreasing energy, so that these activations come 
primarily from the lower energy range of the neutron 
spectra. Typically 70% of the activations come from 
neutron energiesless than 70 keV, except at the 22.2- 
mm radius where there is a significant 14-MeV con- 
tribution. These activations provide a stringent test 
of the transport codes, because the lower energy 
region of the neutron spectrum is populated through 


TABLE 


successive 
downscattering 
interactions. 
Unfor- 
tunately, the (n,~) reaction cross sections for most of 
the target nuclei used are not well known. 
The results of Calculation IV are given in Table 
XII and shown in Figs. 13-17. Only for gold, whose 
activation cross section is thought to be well known, 
are Calculation IV and the experiment in satisfac- 
tory agreement everywhere. Satisfactory agreement 
for scandium, yttrium, and thulium is also obtained 
everywhere except possibly at the 300-mm position, 
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Fig. 13. 
Ratio 
of observed-to-calculated 
4’Sc(n,y)46Sc 
activation 
us radius. 
The error bars are from 
experimental 
uncertainties 
only. 


XIII 


RATIO 
OF OBSERVED-TO-CALCULATED 
‘88U(n,f)”Mo ACTIVATIONS 


Distance From Source 


Calculation 
50.0mm 
76.15mm 
126.0mm 
200.0mm 
300.0mm 


I 
0.880 
0.844 
0.879 
0.925 
1.099 
II 
0.964 
0.908 
0.904 
0.888 
0.988 
m 
1.036 
0.990 
1.000 
0.978 
1.039 
lv 
1.030 
l.ml 
1.004 
0.989 
1.047 
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where there are rather large discrepancies in the 
scandium 
and 
thulium 
results. The 
calculated 
2’8U(n,~) activations are everywhere smaller than 
observed (generally by 10 to 207.), suggesting that 
the activation cross section used is too small, par- 
ticularly in the 14-MeV region, because the dis- 
crepancy is even larger at 22.2 mm. Note that the 
postulated off-center position of the experimental 
source is not expected to influence the (n,~) activa- 
tions at 22.2 mm by more than a few per cent. Final- 
ly, the fact that the observed-to-calculated activa- 
tion ratio is generally large at 300 mm, the outaide 
position, suggests that there may be some room 


return, 
which 
has been 
neglected 
in the calculations. 


A calculation that attempted to simulate the room 
return, however, failed to significantly modify the 
results. Moreover, if sufficient room return were pre- 
sent to explain the scandium and thulium dis- 
crepancies, a large discrepancy would then be in- 
troduced in the gold results at 300 mm because of 
the large value of the Au(n,~) cross section at very 
low energies. 
As for the (n,f) activations, the two Monte Carlo 
calculations are in very close agreement. This is not 
surprising because Calculation III is used to provide 
the lower energy fluences for Calculation IV. The 
two S. calculations are in much poorer agreement 
with the experiment; the more reliable 21-group 
calculation 
gives somewhat 
better results. The 
resultsof the four calculations of the Au(n,~) activa- 
tions, which are typical, are”given in Table XIV. 


D. 
The “’Ir(n,n’)’’’mIr 
Activation 


The cross section for the Issfi(n,n’)igsmbreaction is 
very 
poorly 
known. 
Consequently, 
we consider 
Calculation IV to be in satisfactory agreement with 
the experiment (see Table XII). 


V. CONCLUSIONS 


It seems, on the basis of this experiment, that our 
understanding of neutron transport through OLiDis 
essentially correct. In most cases, discrepancies be- 
tween calculation and experiment can be plausibly 
explained in terms of imprecisely known cross sec- 
tions. The cross-section data relevant to transport 
through ‘LiD 
appear 
to 
be 
essentially correct, 
although there is some indication that the 14-MeV 
inelastic or nonforward elastic cross sections on 
deuterium and/or *Li may be a few per cent too 
small. The cross-section data used for the (n,2n) and 
(n,f) activations seem to be basically correct. Ex- 
cept for mu(n,~), 
the (n,y) activation cross sections 


also seem to be reasonably accurate, even though 
only the cross sections for gold are well known ex- 
perimentally. The ‘Li(n,n’a) 
cross-section data are 
quite possibly in error (too high) by -15% in the 14- 
MeV region, whereas the ‘Li(n,a) 
cross sections 
seem to be quite accurate. 
For neutron-transport computations, the Monte 
Carlo methods provide more realistic resultsthan do 


TABLE 
XIV 


RATIO 
OF OBSERVED-TO-CALCULATED 
‘WAu(n,y)’”Au 
ACTIVATIONS 


Distance From Source 


Calculation 
22.2mm 
50.0mm 
76.15mm 
126.0mm 
200.Omm 
300.0nun 
— 
— 


I 
0.752 
0.740 
0.726 
0.740 
0.803 
0.934 
II 
0.936 
0.926 
0.883 
0.847 
0.848 
0.924 
III 
1.008 
1.067 
0.993 
0.964 
0.938 
1.048 
Iv 
1.026 
1.059 
1.041 
0,983 
0.964 
1.079 
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the S. techniques, although they require signifi- 
cantly more machine time. It is reassuringthat our 
most 
realistic 
simulation 
of 
the 
experiment, 
Calculation IV, provides the best overall fit to the 
data. We believe that further improvement in the 
calculation could be obtained only with the incor- 
poration of improved nuclear data. 
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